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Abstract. Likely astrophysical sources of detectable high-energy (>> TeV) neutrinos are 
considered. Based on gamma-ray emission properties, the most probable sources of neutrinos 
are argued to be GRBs, blazars, microquasars, and supernova remnants. Diffuse neutrino 
sources are also briefly considered. 



1. Introduction 

This discussion of sources of high-energy astrophysical neutrinos focuses on candidate discrete 
sources of neutrinos (see also [lj). Arguments for detectable neutrino sources are based on the 
7/1/ connection: photohadronic and secondary nuclear processes, which are the most important 
astrophysical neutrino production mechanisms, will also produce 7 rays. By identifying the 
brightest and most fluent 7-ray sources, we therefore identify the most likely neutrino point 
sources. This argument is, however, far from airtight. A source can be 7-ray bright without being 
neutrino bright if the 7-rays originate from leptonic processes, for example, when relativistic 
electrons Compton scatter ambient target photons to 7-ray energies. Conversely, a bright 
neutrino source can be 7-ray dim if the 7 rays are attenuated. This can occur either in a 
buried source, where the surrounding material or the local radiation field provides large opacity 
to 7 rays or, for cosmologically distant objects, when high-energy 7-rays are attenuated by the 
extragalactic background light (EBL). 

At EGRET and GLAST energies (E rs 100 MeV - 10 GeV), attenuation by the EBL is 
unimportant, even for the highest redshift objects. Thus the EGRET catalog can be enlisted 
to identify the brightest 7-ray sources and, by the 7-1/ argument, the most probable neutrino 
point sources. The new discoveries with HESS, MAGIC, and VERITAS at TeV energies provide 
additional important information, especially for galactic sources. 

Diffuse 7-ray emissions include those that are genuinely diffuse, namely cosmogenic GZK 
neutrinos and neutrinos produced by cosmic-ray interactions in our Galaxy, as well as apparently 
diffuse neutrino fluxes from the superposition of unresolved sources, such as star-forming galaxies 
and clusters of galaxies. 



2. The 7-z/ connection, 7-ray sources, and candidate neutrino sources 

Photohadronic and secondary nuclear production processes are the two principal mechanisms 
responsible for high-energy neutrino production. The dominant near-threshold channel for 
neutrino production through either of these processes involves the excitation of a A + isobar. 



For proton/photon and proton/nucleon interactions, these reactions take the form 
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where 7' represents a target photon and N represents a target nucleon. The charge-changing 
reaction takes place one-third of the time. For a large number of such reactions, therefore, the 
outcome is a high-energy lepton and three high-energy neutrinos for every four high-energy 7 
rays. In addition, neutron /3-decay forms an electron and a neutrino with ~ 10 2 times less 
energy than the other secondaries. Approximately four times as much energy is radiated in 
electromagnetic secondaries as in neutrinos. Any high-energy neutrino source will therefore be 
a strong 7-ray source if the 7 rays reach the observer without intervening attenuation. 

For such 7-ray sources, the "best-bet" neutrino sources are those where a neutrino telescope 
such as IceCube can detect a neutrino above the cosmic-ray induced background. IceCube's 
background is at the level ^ QAv{> 1 TeV)per yr per square degree and 0.08^(> 10 TeV) per 
yr per square degree [2]. To be detected with a km-scale neutrino telescope such as IceCube, 
the neutrino fluence, and therefore the photon fluence, must be at the level of <I> 7 ^ 10 ergs 
cnr 2 , as we now show [3J. 

The detection probability for muon neutrinos is P u (614) « 10~ 4 ei4 , for 0.1 eu 10 , 
where £14 = e^/100 TeV = e^/160 ergs |4j. Consider a neutrino source with number flux 
4>v{^ia) = dN/dAdtdei^ = Ke^l in the energy range 0.1 ^ £14 ^ 10. The neutrino energy 
fluence received during the time interval At is <3?2,(ergs cm -2 ) = (160 ergs) At J 1( J de\A en ^(£14). 
The detection of N v muon neutrinos occurs when N v = AAtJ 01 deu ff„(ei4)</v(ei4) = 
10~ 4 A <3?j,(ergs cm~ 2 )/160 ergs > 1 , where A is the detector area. Thus the detecttion of 
one v„ with a km-scale neutrino detector (A = 10 10 cm 2 ) requires that the neutrino fluence 
&u ~ 10~ 4 ergs cm -2 . By the 7-^ connection, this means that the 7-ray fluence from a source 
must be <3? 7 S> 10~ 4 ergs cm~ 2 . 

EGRET reported results of their observations in units of W~ s 4>-$ ph(>100 MeV) cm" 2 s _1 
[5]. For a flat vF v spectrum, the mean photon energy in the EGRET range, ~ 100 MeV - 
5 GeV, is ~ 400 MeV, which corresponds to a bolometric vF v flux in this energy range of 
~ 6 x 10~ 12 (/>_g ergs cm -2 s . The integral photon flux sensitivity for a two-week on-axis 
pointing with EGRET is at the level of 0_g ~ 15, so that the limiting sensitivity of EGRET for 
a 5a detection was ~ 10~ 10 ergs cm -2 s _1 . To reach a fluence of 10 _4 $_4 ergs cm -2 therefore 
requires that Atc/>_g = 1.5 x 10 7 $„4 s. The standard EGRET observation lasted for two weeks. 
When taking into account Earth occultation, the typical observing time was w 6 x 10 5 s. Thus 
sources with 0_s 3> 30 fulfill the requirement that if neutrinos are produced with comparable 
fluence as 7 rays, then they would be detectable neutrino sources. Even though the EGRET 
energy range is at a much lower value than that of the neutrino telescopes, the 7 rays produced 
in association with the neutrinos could cascade down into this energy range. Extrapolating the 
7-ray spectrum into the PeV range means that bright 7-ray sources with nearly flat vF v spectral 
indices (or photon number indices a p h ~ 2, defining photon number fluxes 4>-y(e) oc e~ aph ) are 
good candidate neutrino sources. 

By examining the Third EGRET catalog [5j, one finds many 7-ray sources that fit this 
criterion. The following objects have at least one and sometimes many two-week observing 
periods during which 0_ 8 > 100: the blazars PKS 0208-512, PKS 0528+134, NRAO 530, 3C 
279, PKS 1622-297; pulsar 1706-44 and the Crab, Vela and Geminga pulsars; the EGRET 
sources associated with the supernova remnants W44, IC 443, and 7 Cygni; the sources 3EG 
J1824-1514 and 3EG J0241+6103 associated with the microquasars LS 5039 and LSI +61 303, 
respectively; and several unidentified EGRET sources, including some in the Galactic plane and 
some at high latitude. In addition to these persistent or flaring sources, a bright Solar flare 



was detected with EGRET, and several GRBs were strongly detected with the EGRET spark 
chamber. We discuss these various classes of sources in increasing likelihood of being detectable 
neutrino sources. 

2. 1 . Solar Flares 

The June 11, 1991 Solar flare radiated a > 100 MeV energy fluence exceeding 10 -4 ergs cnr 2 [6j. 
The flare spectrum was fit by a slowly decaying (~ 255 minutes) pion emission component and a 
fast-decaying (~ 25 minutes) electron bremsstrahlung component. The very soft 7-ray spectrum, 
with number index £3-4, and lack of evidence for ^> GeV proton and ion acceleration in Solar 
flares (for example, from ground-based neutron monitors) make it unlikely that very high-energy 
neutrinos are produced by Solar flares, though they could be sources of GeV - TeV neutrinos. 

2.2. Pulsars and Pulsar Wind Nebulae 

Pulsars are the brightest point sources in the EGRET catalogs, with > 100 MeV vF v fluxes 
of some pulsars exceeding 10 -9 ergs cm" 2 s -1 . Moreover, they are persistently bright. 
Consequently it only took EGRET ~ 1 day to measure pulsed fluences £ 10 -4 ergs cm -2 . 
It is unlikely, however, that they would be detectable neutino sources. The pulsed spectrum 
always cuts off below several hundred GeV. If this is due to an electromagnetic cascade, as 
expected in both polar cap or outer gap models for pulsed 7-ray emission, then the emission 
originates from electron acceleration and cascades, which would not produce neutrinos. 

The nebulae formed by the cold but highly relativistic MHD winds expelled by rotating, 
highly magnetized neutron stars could accelerate protons and ions that would subsequently 
undergo interactions to produce neutrinos. The Crab nebular emission is at the level of 10 -10 
ergs cm -2 s _1 in the EGRET energy band, but is convincingly explained as the self-Compton 
component from the electrons that radiate the synchrotron nebular emission. 

The pulsar wind nebulae discovered at TeV energies with HESS could be the result of 
accelerated protons that interact with ambient material to form the power-law spectra measured 
from these sources. In its preliminary galactic plane scan [7], HESS reached 5a sensitivities at 
the level of 3 x 10~ n ph(> 100 GeV) cm~ 2 s _1 implying, for a mean photon energy of 400 
GeV, a limiting bolometric sensitivity of 2 x 10~ n ergs cm" 2 s _1 . For example, the pulsar wind 
nebula in the supernova remnant MSH 15-52 radiated ~ 3.3 x 10 -11 ergs cm -2 s _1 in the energy 
range 0.3 - 40 TeV with a p h = 2.27 ± 0.2 [8]. Extrapolating into the PeV range gives a flux 
of ~ 5 x 10 -12 ergs cm -2 s _1 , or a fluences £ 10~ 4 ergs cm" 2 in ~ 1 years time. Because of 
similar morphology to the keV synchrotron emission, the TeV 7 rays are likely to be the result 
of Compton-scattering rather than hadronic emission. If this emission is the result of hadronic 
processes, however, such sources could be marginally detectable neutrino sources (see [SlIIO] for 
a more on Galactic neutrino sources). 

2.3. Supernova Remnants 

Some of the brightest EGRET sources are associated with supernova remnants |llj . and they 
display relatively hard spectra (~ —2). They are also believed to accelerate cosmic rays to the 
knee of the spectrum (~ 3 PeV). If due to 7r° decay emission from secondary nuclear production, 
the spectra will soften at energies well above the ir° peak, which is at several hundred MeV in 
a uF v representation. In the TeV energy range, a bolometric vF v flux reaching 10~ 10 ergs 
cm -2 s _1 is observed with HESS [12] from RX J1713. 7-3946. One would optimistically expect 
[T3] this SNR to be detected as a neutrino sources after only some months of observing except, 
unfortunately, the 7-ray spectrum exhibits a remarkable cutoff above ~ 10 TeV, though IceCube 
should still detect one or two v^{> 1 TeV)/yr [10] . The cutoff might indicate that acceleration 
to higher energies has not yet occurred, though the remnant should be well into the Sedov phase. 
The higher energy cosmic ray protons and ions, being more diffusive, could leave the acceleration 



region quickly. The emission could also be leptonic, arising from Compton-scattered ambient 
radiation. 

2-4- Microquasars and X-ray Binaries 

The discovery that LS 5039 is a TeV source |14| confirmed the association of the EGRET source 
3EG J1824-1514 with LS 5039 [T3]. Its bolometric vF v flux in the 0.2 - 10 TeV range is at 
the level of ~ 10 -11 ergs cm" 2 s _1 . The flux is modulated at the 3.9 day orbital period |16j . 
Orbital modulation was also recently reported by the MAGIC collaboration from the northern 
hemisphere microquasar LSI +61 303 [T7j, with a mean bolometric uF v flux of ~ 4 x 10 -11 ergs 
cm -2 s _1 in the 0.2 - 4 TeV range and mean photon spectral index a p h = 2.5. 

Both of these sources are high-mass microquasars: the companion stars in LS 5039 and LSI 
+61 303 have masses ~ 23M and ~ 10M & , respectively. The orbital modulation shows that 
the TeV emission has to be produced in the vicinity of the binary system. Both leptonic [181 [T9] 
and hadronic models [20J for microquasar jet emission have been proposed. Even though their 
spectra are soft or cut off at multi-TeV energies, this could be a result of 77 attenuation, so that 
the actual 7-ray and neutrino production spectrum from hadronic interactions could extend to 
very high energies [20|, [2T] . 

In addition to these two microquasars, the binary system B1259-63, consisting of a pulsar 
and a high-mass Be star, is a TeV source [22] , This raises the interesting possibility [23J that 
all three sources are binary systems consisting of a pulsar and a high-mass star, with the TeV 
emission due to nonthermal particles accelerated by the shock formed by interactions between 
the MHD wind of the pulsar and the stellar wind of the high-mass star. The formation of the 
emission in the inner region of the system would naturally follow from this scenario. Accreting 
X-ray binaries have also been proposed [23] as detectable neutrino sources when protons and 
ions, accelerated in the magnetosphere of the system, collide with material of the accretion disk, 
though the absence of TeV emission suggests they would not be bright neutrino sources. 

2.5. Blazars 

The intense, highly variable 7-ray fluxes from blazars suggest that they are also bright neutrino 
sources. It is important to distinguish between the flat spectrum radio quasars (FSRQs), which 
have strong atomic emission lines in their spectra, from the BL Lac objects with their nearly 
featureless continua. All known TeV blazars are X-ray selected BL Lac objects. 

Though one might think that TeV blazars are the most probable neutrino sources, given 
that particle acceleration to 3> TeV energies must take place in these sources, it is more likely 
that FSRQ blazars are bright neutrino sources [25]. First, the brightest FSRQs have vF v fluxes 
~ one order of magnitude brighter than the TeV vF v fluxes of the brightest BL Lac objects. 
The absence of TeV radiation in FSRQs is probably a consequence of 77 attenuation on the 
EBL, and does not indicate that TeV 7 rays are not produced. Moreover, the intense scattered 
accretion-disk radiation in the vicinity of the supermassive black hole in FSRQs, as indicated 
by the strong atomic lines, provides an important source of target photons for photohadronic 
production [26]. Making the conservative assumption that the energy injected in protons is 
equal to the energy inferred from observations of the electron synchrotron radio/X-ray emission, 
we [251 l2"B] have shown that IceCube could detect one or several neutrinos during bright FSRQ 
blazar flares, such as that observed from 3C 279 in 1996. 

2.6. Gamma Ray Bursts 

This topic was recently reviewed [3], and details can be found there. Neutrino production in 
GRBs depends most sensitively on two parameters: the baryon-loading factor and Doppler factor 
of the GRB blast wave. The baryon-loading factor refers to the energy in nonthermal protons 
compared to the electromagnetic energy inferred from direct measurements of the keV/MeV 
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Figure 1. GZK neutrino flux using a star 
formation rate history used to model the 
UHECR spectrum from GRBs [35]. 
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Figure 2. GZK neutrino flux using the star 
formation history used to fit the statistics of 
GRBs with measured redshifts l36l. 



emission from GRBs. Provided that the baryon-loading factor is 3> 10, which is required if 
GRBs are the sources of UHECRs, and the Doppler factor is ^ 200, neutrinos from GRBs are 
detectable with IceCube or a northern hemisphere neutrino detector. These calculations [29] are 
made in the framework of the collapsar model, with values of the Doppler factor in the range 
commonly expected for GRB outflows. Anomalous 7-ray emission components |27} 128] in GRBs 
give further evidence for hadronic acceleration by GRB blast waves. 

3. Diffuse Neutrinos 

Cosmic ray interactions in the disk of the Milky Way will make a diffuse neutrino glow |30j . 
A flux of ps 160 v (> 250 TeV) km" 2 yr" 1 from a 5 square degree region surrounding the 
galactic center is calculated in Ref. PQ, though the exact result is sensitive to the hardness of 
the cosmic-ray spectrum. Superpositions of emissions from the classes of point sources listed 
above, in particular, GRBs and blazars, will make a diffuse high-energy neutrino background 
radiation. Dim extragalactic 7-ray sources can also make important neutrino backgrounds due 
to their abundance. 

Star-forming galaxies, which include normal spiral galaxies, starburst galaxies, and infrared 
luminous galaxies, might be considered as likely neutrino point sources because cosmic rays 
would certainly be accelerated by the core-collapse supernovae resulting from the late stages 
of evolution of the massive stars in these systems. These systems are, however, relatively dim 
7-ray sources and, furthermore, have soft spectra. The only extragalactic galaxy that was 
detected with EGRET was the Large Magellanic Cloud [31], with an integral 7-ray flux equal 
to 14.4(±4.7) x 10" 8 ph(> 100 MeV cm s and a spectral index of s = 2.2, implying a vF v 
flux of ps 2.3 x 10- n (£yi00 MeV) -0 - 2 ergs cm" 2 s" 1 . Thus it would take > 2 years to detect 
neutrinos from the LMC unless there was an anomalous hardening of the spectrum. 

Nevertheless, the superpositions of the neutrino emissions from star-forming galaxies will 
form a guaranteed background. An estimate [35] of the neutrino background based on the 
synchrotron radio luminosity associated with cosmic-ray acceleration in star-forming galaxies 
is at a level detectable by IceCube, though the assumptions and derived intensity have been 
challenged [55] . 

Photohadronic interactions of UHECRs with photons of the EBL forms the diffuse cosmogenic 
GZK neutrino flux. The spectrum of this background depends on the UHECR activity in the 
early universe |34| . Figs. 1 and 2 show calculations (work in preparation with J. Holmes) of the 
effect on the GZK neutrino background due to different histories of GRB production. 



4. Summary 

The best-bet neutrino sources are GRBs, blazars, microquasars, and SNRs because of their bright 
and hard 7-ray spectra that could originate from hadronic processes. These results should be 
placed within the context of theories of cosmic ray origin. Supernova remnants are thought to 
be the sources of cosmic rays to at least the knee of the cosmic ray spectrum, because they have 
adequate power and will produce strong shocks. But the TeV spectrum of RX J1713. 7-3946 
shows a cutoff at ~ 10 TeV. Supernova remnants probably differ greatly in their cosmic ray 
acceleration efficiencies, with Type lb/c supernovae and those associated with GRBs being the 
strongest such accelerators. Microquasars and Be/X-ray binaries, though less powerful than 
SNRs, will also form strong shocks and will have stellar wind material that can provide a target 
for secondary nuclear production. 

UHECRs must be extragalactic given the ~ fiG strength of the Milky Way's magnetic field, so 
that GRBs and blazars, the brightest extragalactic 7-ray sources, are the most probable sources 
of UHECRs and high-energy neutrinos. The highly relativistic outflows and shocks required to 
model these systems can accelerate particles to « 10 20 eV. 
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